Soliton transverse instabilities in nonlocal nonlinear media 
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We analyze the transverse instabilities of spatial bright solitons in nonlocal nonlinear media, both analytically 
and numerically. We demonstrate that the nonlocal nonlinear response leads to a dramatic suppression of the 
transverse instability of the soliton stripes, and we derive the asymptotic expressions for the instability growth 
rate in both short- and long-wave approximations. © 2008 Optical Society of America 
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Symmetry-breaking instabilities have been studied in 
different areas of physics, since they provide a simple 
means to observe the manifestation of strongly non- 
linear effects in nature. One example is the transverse 
instabilities of spatial optical solitons [1] associated 
with the growth of transverse modulations of quasi-one- 
dimensional bright and dark soliton stripes in both fo- 
cusing [2-4] and defocusing [5] nonlinear media. In par- 
ticular, this kind of symmetry- breaking instability turns 
a bright-soliton stripe into a array of two-dimensional 
filaments [6] and bends a dark-soliton stripe creating 
pairs of optical vortices of opposite polarities [7]. Con- 
sequently, the transverse instabilities set severe limits 
on the observation of one-dimensional spatial solitons in 
bulk media [8]. 

Several different physical mechanisms for suppressing 
the soliton transverse instabilities have been proposed 
and studied, including the effect of partial incoherence 
of light [9, 10] and anisotropic nonlinear response [10] in 
photorefractive crystals, and the stabilizing action of the 
nonlinear coupling between the different modes or polar- 
izations [11]. Recently initiated theoretical and experi- 
mental studies of nonlocal nonlincaritics revealed many 
novel features in the propagation of spatial solitons in- 
cluding the suppression of their modulational [12] and 
azimuthal [13] instabilities. In this Letter we demon- 
strate that significant suppression of the soliton trans- 
verse instabilities can be achieved in nonlocal nonlinear 
media, and we derive analytical results for the instabil- 
ity growth rate in both long- and short-scale asymptotic 
limits. 

We consider the propagation of an optical beam in a 
nonlocal nonlinear medium described by the normalized 
two-dimensional nonlinear Schodinger (NLS) equation, 
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optical refractive index, and the parameter d stands for 
the strength of nonlocality. The model ([I]) describes the 
light propagation in different types of nonlocal nonlinear 
media including nematic liquid crystals [14]. 

We look for stationary solutions of Eq. (p} in the form 
of the bright soliton stripes, E{x,y;z) = u(x) exp(i/3z), 
where u(x) is a (numerically found) localized function, 
u(±oo) = 0, and (3 is the (real) propagation constant. 

The transverse instability of quasi-one-dimensional 
solitons in nonlocal nonlinear media is investigated by 
a standard linear stability analysis [1], by introducing 
the perturbed solution in the form: n = no{x) + eSn, and 
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where (u , n ) is the solution of Eqs. (JTJ), e 1 is a small 
perturbation, and v(x), w(x), and 5n(x) are perturbed 
amplitudes that are modulated in the transverse y direc- 
tion with the wavenumber p. The instability growth rate 
is defined as an imaginary part of the eigenvalue A. 

Substituting these asymptotic expansions into Eq. (fT]), 
in the first order of e we obtain a set of linear equations, 
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Sn = d ( - p 2 ) Sn + 2u v, 



where Aj_ = d 2 /dx 2 + d 2 /dy 2 , E = E(x,y;z) is the 
slowly varying electric field envelope, n = n(x, y; z) is the 



which we then study numerically and analytically. 

Typical solutions for bright solitons of the model ([1]) 
are shown in Fig.QJa) for local (solid curve) and nonlocal 
(dashed and dashed-dotted curves) media. Figure [2(b) 
shows the growth rate of the soliton transverse instabil- 
ities vs. the modulation wavenumber p, for local [d = 0) 
and nonlocal {d = 0.5 and d — 1) nonlincaritics, re- 
spectively. It can be seen that nonlocality reduces the 
growth rate of the transverse instability of the soliton 
stripe. Moreover, the cutoff transverse wavenumber p c 



of the gain spectrum becomes smaller as the value of the 
nonlocality parameter d grows. 

To describe the suppression of the soliton transverse 
instabilities quantitatively, we calculate the dependence 
of the cutoff transverse wavenumber p c and the maxi- 
mum growth rate on the strength of nonlocality. d, shown 
in Figs. [1Jc,d). We observe that the maximum growth 
rate decreases significantly at large values of the nonlo- 
cality parameter, and eventually it approaches zero when 
d — > oo. The cutoff wavenumber p c of the transverse in- 
stability domain becomes smaller as the value of nonlo- 
cality grows. In the limit of very large values of d, the 
cutoff wavenumber vanishes as well. Consequently, the 
soliton stripes tend to become more stable when the de- 
gree of nonlocality increases. 

Next, we analyze the transverse instability of bright 
solitons in nonlocal nonlinear media by applying a vari- 
ation method, in accord with the following steps. First, 
we expand the nonlocal refractive index function into se- 
ries in the nonlocality parameter [15] in the terms involv- 
ing the refractive index in Eqs. @. Second, we apply the 
method of Ref. 3 to construct the asymptotic expansions 
of the elliptical problem defined by Eqs. ((21) ■ We employ 
the corresponding ansatz defined as: v = Vq + Tv\ and 
w = wq + Twi, where T is an imaginary part of the 
eigenvalue of the linear stability problem. 

In the long-scale expansion, we use the following so- 
lutions: vq = 0, vi = rf 2 [sech(6ia;)]/(ia; 2 , and wq = uq, 
where the parameter b\ is obtained by minimizing the 
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Fig. 2. Growth rate of the soliton transverse instabilities 
defined from numerical simulations and variational ap- 
proach. Solid and dashed lines are numerical data for lo- 
cal (d = 0) and nonlocal {d = 0.5) nonlincarities. Upper 
and lower triangular marks: long-scale and short-scale 
asymptotic expansions for local nonlinearity. Squares 
and circles: long-scale asymptotic expansion based on 
the linearized nonlocal eigenvalue system expanded up 
to the first and second orders in d. Dots and crosses: 
the short-scale asymptotic expansion up to the first and 
second orders in d. 
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Fig. 1. (a) Intensity profiles of the bright solitons in 
local (d = 0, solid) and nonlocal [d — 0.5, dashed; 
d = 1, dashed-dotted) nonlinear media, (b) The insta- 
bility growth rate for bright solitons vs. the transverse 
wavenumber p for local (solid) and nonlocal (d = 0.5, 
dashed; d = 1, dash-dotted) nonlinear media. (c,d) Cut- 
off value of the transverse wavenumber and the maxi- 
mum growth rate vs. the nonlocality parameter d. 



system Lagrangian corresponding to Eq. |[2J , 
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Wavenumber j3 and solution uq = Asech(aa;) are 
obtained by minimizing the Lagrangian of the non- 
local NLS equation based on the same expansion 
technique. For a given power P and nonlocality pa- 
rameter d, we obtain (3 = (l/6)Pa + (l/6)a 2 + 
(2/15)dPa 3 , A = (Pa/2) 1 ' 2 , and a = [-5 + (25 + 
60dP 2 ) 1/2 ](12dP)" 1 . While evaluating the system La- 
grangian, we use the functional expansion, sec(ax) « 
sech(6i.T) — sech(6ix)tanh(&ia;)(a — b\)x to deduce a close 
form for the Lagrangian in Eq. ([3|). The parameter b\ 
can therefore be obtained as a function of both d and P. 
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Fig. 3. Evolution of a modulated bright soliton stripe 
in (a-c) local and (d-f) nonlocal nonlinear media at the 
distances z = 1.0, 6.0 and 10.0, respectively. 

Then, the instability growth rate is obtained as 

1 



r- 



(w ,u An(vi;p = 0)) 2 
— V 



(wo,wi) 



6n(vi;p = 0) 



1 + dp 2 



(1 + dp 2 ) 2 



(4) 



where we define (u>o,Vi) = J_ 6xWqVi. 

For the asymptotic expansion near the cutoff 
wavenumber p c , we use the ansatz vq — scch(6o^c), 
wo = 0, and w\ = sech(cix). Again, to obtain bo we min- 
imize the Lagrangian L Vl as a function of the wavenum- 
ber p. Then the cutoff wavenumber p c is found from the 
condition L Vl = with an optimized value of 6o- The 
resulting parameters are deemed as constant in the La- 
grangian for wi, 
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media, the long-scale expansion to the first order in d 
(circles) is insufficient, but the second-order expansion 
(squares) gives a good result. In the short-scale expan- 
sion, both the first- (dots) and second-order (crosses) 
expansions arc in a good agreement with the numerics. 

To observe the consequence of the predicted instability 
suppression, we study numerically the evolution of the 
soliton stripe described by Eqs. ([1]). Figures [3]^ a-f) com- 
pare the transverse instability of a bright-soliton stripe 
with the power density 2 in local and nonlocal media, 
with the initial field modulated transversally with the 
maximum growth rate. In Fig. [3ja-c) we show the snap- 
shots of the soliton stripe evolution in a local nonlinear 
medium at z = 1.0, 6.0 and 10.0, respectively. We ob- 
serve that at z = 10.0 [Fig.[3](c)], the bright-soliton stripe 
decays into a sequence of filaments due to the modulation 
in the y direction. In a sharp contrast, there is no visible 
decay of the soliton stripe in a nonlocal medium (d = 1) 
[Fig. [3j d-f)], therefore confirming our major conclusion 
that the soliton transverse instabilities are suppressed 
substantially in nonlocal nonlinear media. 

In conclusion, we have analyzed the transverse insta- 
bilities of spatial solitons in nonlocal nonlinear media. 
By employing the linear stability analysis and numerical 
simulations, we have demonstrated that nonlocal nonlin- 
ear response can suppress significantly the transverse in- 
stabilities allowing experimental observations of the sta- 
ble propagation of the soliton stripes in nonlocal media. 
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to deduce parameter c\ as a function of P and d. Simi- 
larly, in the short-scale expansion the amplitude growth 
rate can be calculated as 

r2 _ (vo,u An(w 1 ;p = 0)) 2 _1 4 

P A P 1 



Sn(wi;p = 0) 



1 + dp 2 



(1 + dp 2 ) 2 



(6) 



Figure [2] compares the analytical results with the re- 
sults of our direct numerical calculations. For nonlocal 
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